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1 Introduction 



The effective action T[0], the generating functional of the one-particle irre- 
ducible (1PI) vertex functions, is a very useful tool in quantum field theories. 
It is widely used in the analysis of their vacuum structure [1] (for more re- 
cent reviews see e.g. [2,3]), and its symmetries are often exploited in order to 
establish their renormalization properties [4] . For vanishing external sources, 
it satisfies the equation 



sometimes referred to as the quantum equation of motion, the quantum coun- 
terpart of the classical equation || = 0. 

In fact, at the lowest order of the semi classical expansion, T[0] coincides 
with the classical action £[</>], and this property, when combined with Eq. 
(JTJ), naively suggests that T[(j)} should be regarded as the quantum action 
of the system, and Eq.(JTJ) as the corresponding quantum dynamical equa- 
tion. This interpretation, however, is correct only under certain conditions. 
Moreover, even when this interpretation applies, Eq. ((TJ) in general cannot be 
straightforwardly related to any semiclassical expansion. 

The dynamical content of Eq.Q was already briefly discussed in [5], 
where a solid theoretical background to the work presented in [6] was pro- 
vided along with the correction of some mistakes and misunderstandings. 
One of the purposes of the present paper is to present a more detailed de- 
scription of the approach followed in [5], with the aim of clarifying the limits 
of applicability of this procedure. 

As it is well known, an ansatz for T[(f>] which is well suited for the ap- 
plications is given by the gradient expansion, an expansion in powers of the 
field derivatives : 



The lowest order approximation, the so called Local Potential Approx- 
imation (LPA), is obtained from Eq.(J2J) once we neglect the higher order 




(2) 
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derivatives and keep the wave function renormalization factor constant, i.e. 
field independent (without loss of generality we can choose Z = 1) : 



m = J d 4 x(-K//(0(*)) + \d^{x)d^{x)) • (3) 
Within the LPA, Eq.((TJ) becomes: 

ouo , 9V ef f((f>) . . 

^ = 30 ' (4) 

i.e. it takes the same form of the classical equation of motion where the 
classical potential is replaced by the effective potential. 

This equation has been widely used in the past as the quantum equation of 
motion to describe the time evolution of the scalar field expectation value [7]. 
More recently, however, some of the limitations in the use of Eq. (@J), and 
more generally of Eq. as a dynamical equation have been noted [8] and 
a different formalism to deal with dynamical problems has been developed 
[8-10]. 

In the present paper we complement these previous studies and push the 
investigation of Eq. (@J) a step further. We carefully analyse the limitations 
in the application of this equation to dynamical problems and show that 
there are physically relevant cases where it can be appropriately used to 
describe the evolution of expectation values. It will turn out that two points, 
overlooked in the past, play a crucial role: the correct consideration of the 
boundary conditions encoded in the definition of the effective action and the 
use of a non perturbative approximation for the effective potential in Eq. (@J). 

Let us proceed now to our systematic analysis of the subtle points related 
to the use of Eq. (Q) as a dynamical equation. As it will be clear in a moment, 
it is essential to specify: 

(i) the framework in which Eq. (0) is derived, to avoid confusion on the 
physical meaning of (p; 

(ii) the choice of the boundary conditions (consistent with point (i)) associ- 
ated to the differential equation (JTJ); 

(iii) the approximation in which T[(p] is computed, that has obviously to be 
well suited for the particular problem considered, in order to get physically 
meaningful results. 

Concerning point (i), one has to be careful in considering cf) as the field ex- 
pectation value because in general this is not true. The effective action T[cf>] 
is the Legendre transform of W[«7], the generator of the connected Green 
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functions, defined by: e^ w[J] = Z[J] = (0,t = +oo|0,t = -oo)j. In the 
Schrodinger picture, the vacuum persistence amplitude, Z[J], can be explic- 
itly written as: 

Z[J\ = (0|t>(+oo,-oo)|0) = (0|T(e~*J-«* (d ~ J(t) * W) )|0), (5) 

where U(+oo, — oo) is the time evolution operator in the presence of the 
source J, and |0) is the ground state of H, the Hamiltonian of the system. 
From Eq. (0), and from the definition of the classical field, (f>(x,t) — sw ^ 
we immediately see that: 

(~,t\ + ,t) 



6J{x,t) ■ 



»(g,«) = < -, t| T |H ;- t> , (6) 



where 

(-,t\ = (0\U(+oo,t) and |+, t) = U(t, -oo)|0) . (7) 

Clearly \+,t) and | — ,t) are (in general) different states. As a consequence, 
4>(x, t) is not (in general) a diagonal matrix element of the quantum field $(x). 
Moreover it may be complex-valued. This is not surprising as T generates 
the 1PI vertex functions that are in general complex quantities satisfying 
causal (Feynman) boundary conditions. Only for particular choices of the 
external source J can we have \+,t) = \—,t); under these conditions, Eq. (JTJ) 
describes the dynamical evolution of the field expectation value of $. 

Effective actions defined according to Eq. (jSj) are called Schwinger-De 
Witt (or "in-out") effective actions [11]. They incorporate boundary condi- 
tions that are appropriate to problems where the transition from asymptotic 
states in the past to asymptotic states in the future is considered (as is the 
case of a scattering process). In passing we note that more general defini- 
tions of "in-out" effective actions can be obtained if the |0, t = — oo) and 
|0, t = +oo) asymptotic states are replaced by more general |in) and |out) 
states, i.e. by more general asymptotic boundary conditions. 

A functional formalism that is in general appropriate to dynamical prob- 
lems is the so called "in-in" or closed-time-path formalism [12-14], where one 
can construct "in-in" effective actions that generate the dynamical evolution 
of expectation values. This formalism is the one most largely used nowdays 
in the applications to inflationary cosmology. 
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In this paper, however, we shall not deal with this latter approach and 
its most recent developments. Our aim is rather to show that, under certain 
conditions, we can still define an equation of motion in the in-out formalism. 

In view of the fact that Eq. (JTJ) is sometimes used as the quantum coun- 
terpart of the classical equation of motion, this problem is certainly worth 
to study. 

Even in those cases where this formalism can be properly employed, how- 
ever, we still have to face the problem mentioned in point (ii). Equation 
(J5J) contains two asymptotic conditions at the initial (t = — oo) and final 
(t = +00) times. These conditions are encoded in the definition of the ef- 
fective action itself, as well as in the definition of the classical field (the 
argument of T). It is then clear that we cannot freely choose certain bound- 
ary conditions for the field and its derivatives at a given initial time t = to, 
together with a "physically convenient" form for the (functional) wave packet 
of the system at the same time, and then evolve the expectation value of the 
quantum field according to Eq. (0). In fact arbitrary initial conditions for (p 
and for the wave packet in general are not compatible with the asymptotic 
conditions that enter the definition of V. 

In the following, we provide arguments showing that, by considering a 
specific set of initial conditions for the expectation value of the field and of 
its derivatives, it is possible to find the form of the initial (i.e. at t = t ) wave 
packet that is compatible with the asymptotic conditions encoded in the def- 
inition of the effective action, and whose dynamical evolution is governed by 
Eq. Our arguments will be strongly supported by the numerical results 
that we shall obtain for the dynamical evolution of the position operator 
expectation value in the quantum mechanical double-well potential problem. 
In this case we are also able to consider the exact dynamical evolution with 
the help of the time-dependent Schrodinger equation, so that we can compare 
the results obtained with Eq. (0) with the exact ones. As we shall see, when 
the initial wave packet is chosen according to our criterion, we find excellent 
agreement between the two results. 

We have found it convenient to investigate the question of the determi- 
nation of the initial wave packet within the framework of the variational 
definition of the effective action [15], that is the generalization to the time- 
dependent case of the well-known variational definition of the effective poten- 
tial [16]. The effective action T[(p] is obtained through the extremization of a 
certain functional of two (a priori) different states, called \ip+,t) and \ijj_,t), 
under the constraint that (tp + , t\&(x)\i/}-, t) = <p(x,t) (see next section). 
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Finally we come to point (iii) and to the problem of the approximations 
employed to compute the effective action and effective potential. The most 
simple approximation of T is given by the LPA, i.e. by Eq.Q, and the most 
straightforward approximation of the effective potential V e ff is given by the 
one-loop potential Vu. This last approximation however presents a serious 
drawback. The exact effective potential (actually the effective action) is a 
convex function of its argument [2,17,18]. However, when the classical (bare) 
potential is not convex (these are the physically most interesting cases), 
at any finite order of the loop expansion the approximated effective poten- 
tial does not enjoy this fundamental property. Alternative non-perturbative 
methods of computing the effective action and potential, though, such as 
lattice simulations [19], variational approaches [20], or suitable averages of 
the perturbative results [21], provide the proper convex shape. In addition 
to the methods quoted above, a non-perturbative convex approximation to 
the effective potential, Vrg, is found within the framework of the wilsonian 
renormalization group (RG) equation [22-26]. We have computed the effec- 
tive potential with the help of this RG equation and then inserted Vrg in 
Eq. (ji). The comparison of the results obtained with Eq. (JTJ with those 
obtained with the help of the time-dependent Schrodinger equation shows 
that Vrg, because of its non-perturbative features, provides an excellent ap- 
proximation for the correct "quantum potential" to be used in the "quantum 
equation of motion". For completeness we will also check the inadequacy of 
V\i to describe non-perturbative regimes with our equations of motion. 

The plan of the paper is as follows. In Section 2 we briefly illustrate 
the basic formalism employed in the following. In Section 3 the central 
argument concerning the validity of Eq. ([I]) is discussed, while its application 
to the harmonic oscillator and to the double-well potential are respectively 
addressed in Section 4 and 5. The former example is treated analytically 
whereas the latter is solved by a numerical analysis. Section 6 contains the 
summary and outlook. 

2 Variational definitions of T[0] and V e ff{4>) 

To set up the tools of the following analysis, we briefly review in this section 
the variational principles that lead to the definitions of the effective action 
and the effective potential, referring to [15] and [16] for details. According 
to [15] the effective action is the stationary, time-integrated matrix element 
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of idt-H 

T[(j>]= dt(^,t\(idt-H)\ij; + ,t), (8) 

J —oo 

where the right-hand side is stationary when the two time-dependent states 
\ip±,t) are varied arbitrarily and independently, but with the two constraints 

<^_,£|l>(f)|^ + ,t> = </>(£, t) (9) 

and 

(^_,#+,t> = l, (10) 
and with the asymptotic boundary conditions 

Km fok,*) = |0>, (11) 

t — >=poo 

where |0) is the ground state of H. 

This variation, together with the constraints © and (|10p. is translated 
into the equations: 

+ y d 3 f J(f,t)$(x)J = w(t)\i/)+,t) (12) 

(id t -H + J d 3 x J(x,t)$(x))\^,t) =w*(t)\ip_,t), (13) 

where J(x,t) and u>(t) are the Lagrange multipliers that implement the two 
constraints Q and (JTUJ). The relation between the couple of states \if)±,t) 
and |±,t), introduced in Eqs. (J7J), is shown in [15]. They are related by a 
phase factor given by the time integral of the Lagrange multiplier w(t). 

If we limit ourselves to considering constant (in space and time) field 
configurations, the effective action is reduced (see Eq. (jSJ)) to the effective 
potential, the generator of the 1PI graphs with vanishing external momenta: 

mU= const =-Veff(<P)Jd 4 X. (14) 

As is well known, V e ff(4>) can be obtained by minimizing the expecta- 
tion value of the Hamiltonian among the normalized time-independent states 
which have a field expectation value equal to (p [16], i.e.: 

K//(0)=min^|W>, (15) 
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with the states \ip) subject to the constraints: 



(16) 



and 



<</#> = 1- 



(17) 



The constrained minimum condition in Eq. (JT5|l generates the time-independent 
Schrodinger equation for a modified Hamiltonian: 



where again J and E are the Lagrange multipliers associated to the two 
constraints (fTp^l and (fTTj) . As is clear from Eq. (JT3J) , the ground state(s) of 
H, which is (are) obtained by solving Eq. (fTSj) with J = 0, is (are) associated 
to the minimum (minima) of V e ff(<f)). 

Except for very few cases, the exact form of the effective action cannot be 
determined, and we have to resort to some approximation. A typical ansatz 
for T, very appropriate for our following considerations, is the gradient expan- 
sion, an expansion in terms of the field derivatives. Suitable approximations 
are obtained by considering truncations to a given order. Although in the 
applications we shall limit ourselves to considering the LPA (see sections 4 
and 5), which is the lowest order in this expansion and is the approximation 
typically adopted in dynamical problems, the considerations that follow do 
actually apply to any order of the expansion (see next section). 

3 The quantum equation of motion 

We have already noted that the classical field <p(x, t) (the argument of T) 
is not in general a diagonal matrix element, i.e. it is not the expectation 
value of the quantum field $ in a given state. However, as we shall see at 
the end of this section, under certain conditions the two states \if)+,t) and 
\ip-,t) of Eq. (jSJ) coincide: \ip+,t) = \ip-,t) = \tp,t), in which case 4>(x,t) is 
the expectation value of $(x) in \ip,t). For the moment we assume that we 
are under these conditions. 

Let us now consider the differential equation for <p(x,t), Eq. (JTJ), once 
the effective action is approximated with a derivative expansion truncated 
at the order m, where 2m indicates the highest number of field derivatives 




(18) 
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(in the approximation considered in Eq. (J3J), it is m = 1). In order to get a 
unique solution of Eq. (JTJ) in the time interval [t ,t\\ (where t\ is arbitrarily 
chosen), we need 2m boundary conditions, which we can fix for instance on 
the manifold t = to, with to chosen as the initial time for the evolution of (j). 

To make contact with the variational principle discussed in Section 2, 
we have to consider the evolution of <p(x, t) from t = — oo to t — +oo, and 
therefore we must provide 4>{x,t) also in the time intervals ] — oo,£ ] and 
[£i,+oo[. We can choose the function <j)(x,t) in these intervals arbitrarily, 
provided the asymptotic conditions in Eq. and the proper matching of 
this function with the unique solution of Eq. (JTJ) in the interval [to, h] at the 
times t and ti are taken into account. Once the function 0(x, t) is assigned 
in the whole range ] — oo, +oo[, the constraint on the right-hand side of Eq. 
is defined and we are therefore able to implement the variational principle 
and to determine the corresponding source, which we indicate with J(x,t), 
as well as the state \ip+,t) = \ifj-,t) = \ip,t). 

We note in passing that the source J(x, t) can also be obtained from its 
well-known relation with the functional derivative of the effective action 

™ J(x,t). (19) 



x, t) 

Equation (fT^j) is the generalization of Eq. ((TJ), when a source, linearly coupled 
to the field, is turned on. Obviously the source J(x, t) associated to the field 
4>(x,t) considered above must vanish in the range [toi^i]- 

So, on the basis of the variational principle, we conclude that there is only 
one state \if},t) that is associated to the specific function <fi(x,t) introduced 
above, and in particular that at t — to the boundary conditions considered 
for the time evolution of the field are in fact related to \ip,to). This is an 
important point because, as mentioned in the Introduction, if the constraints 
on the variational principle (see Eqs. 0, (fTUj) and (JUJ)) are neglected, one 
can in general find more than one normalized state that is compatible with 
the set of boundary conditions fixed at t . However, among these states, only 
\i/},t) is truly related to the effective action and therefore to the full evolution 
of 4>(x, t) in the time interval ] — oo, +oo[. 

We can now understand what we have already anticipated in the Intro- 
duction concerning the applications. In a typical cosmological application 
the effective action is truncated as in Eq.fjSJ), a convenient form of the wave 
packet at the initial time t = t (typically gaussian) is taken, and the initial 
values for the expectation value of the scalar field (p an d of its first time 
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derivative are usually taken to be </> = constant and = 0. On the basis 
of the previous considerations, we easily understand that the time evolution 
of (j) is not always related to the time evolution of the wave packet that has 
been considered, or, in other words, <p(x, t) need not be the mean value of 
$(£) in \ip, t). 

Our main problem is then the determination of to), so that we can have 
a direct physical interpretation of the function (p(x, t) as the time evolution 
of the expectation value of the field for that particular state. Even though it 
is certainly not easy to find a solution to our problem for generic boundary 
conditions, we expect that under certain circumstances it should be possible 
to determine \if),to). In the following we consider which is relevant 

to the physical applications (see above), where \ifi,t ) can be found. Namely 
we take at t = t a constant field 4>(x , t ) = (f) with vanishing derivatives 
up to the (2m) th order. Our goal is now to identify the state \i/),to) defined 
through the variational principle. 

To this end we consider a different physical problem, namely the case in 
which, for t > t , the field is a constant : (j) = O - F° r this new problem 
and for — oo < t < t , the source J has to be equal to the one of the previous 
case (J = J(x,t)), while for t > t , J is given by the constant J ? which 
corresponds to O - This constant is obviously fixed by the condition: 

sr[<f>) 



= Jo- (20) 

0=00 

With this assignment for the source J, the solution of Eq. (J12)) for t > t 
must be the time-independent state 

|V(Jo)> = |V^o>- (21) 
In fact in this case Eq. (|T2*|) becomes: 

(H-Jof d 3 f $(f) - £ ) mJo)) = (22) 

where w(t) = —Eq is now a time-independent constant, and the two con- 
straints in Eqs. Q and (fTUj) are: 

(^J )|£(£)|V>(Jo)> = 0o, (HJoMW) = 1- (23) 

Clearly this is nothing else than the time-independent variational prin- 
ciple that defines the effective potential (see Eqs. (f*Pf|) - (frS|0 . We then have 
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found that the state \i[),to) = l^(^o)) we were looking for is nothing but 
the normalized time-independent state that provides the minimum expec- 
tation value of the Hamiltonian H among the normalized states with field 
expectation value equal to 4>o- 

There is one main assumption behind this result, and it concerns the 
source J. In fact, for our purposes, we have considered two different sources, 
which are both equal to J(x,t) in the range — oo < t < t , whereas for 
t > to one is vanishing while the other is a non- vanishing constant J = Jq. 
Note that, by construction, J = J(x,t Q ). Therefore in the former case 
the source is discontinuous at t = to and we are assuming that this has 
no consequences on the determination of the solution of the time-dependent 
Schrodinger equation. This discontinuity could be replaced with a sharp but 
still continuous change of the function J around t — to, but again we would 
have to assume that this has no consequences on the determination of the 
state. 

Finally, let us come back to the important question, mentioned at the 
beginning of this section, that at least in some cases the effective action (or 
an approximation of it) can be obtained as a diagonal matrix element, i.e. 
in Eq. (JHJ) we have |^>_, t) = \ijj + ,t). We now argue that this should be the 
case when the motion of the field 4> in the time interval [to>^i]> i- e - the time 
interval where the motion of cp is governed by Eq. (JIJ) plus a set of boundary 
conditions, is periodic. Suppose that this time interval covers n periods (let 
us call At the period), i.e. that t\ = t + n At. Then, at time t = ti, the 
field 4> and its derivatives take the same values as at t = t . 

Now we can obviously make a shift on the time axis without changing 
the physics and take the point t = to be symmetric with respect to to and 
tx, i.e. —to — t\ — t in the new frame. Let us now choose J in Eq. (J12j) 
to be symmetric for time reflection, namely J must be taken in such a way 
that the ground state at t = — oo is driven onto \ip, —t) at t = — t (and as 
before we indicate this particular source as J(x, t) with — oo <£<—£), then 
J = within the range [— t, t\ (in this interval cf> covers n periods), and, for 
t < t < +oo, the source is turned on symmetrically, i.e.: J(x,t) = J(x, —t). 

Clearly, with such a choice, the evolution of the two states in Eqs. (II 2j) 
and (j!3|) is symmetric for time reversal and therefore equal, i.e. \ip+,t) = 
\ip-,t). It must be noticed that the class of periodic motions of (f> obviously 
includes the case of constant <fi (i.e. the case associated to a constant non- 
vanishing source), which was also considered above. 

The conclusion of this section is that, at least for those cases in which 
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we consider a truncated derivative expansion of the effective action which 
gives rise to a periodic dynamical evolution of 0, this latter quantity has the 
meaning of the expectation value of the field operator, and the associated 
quantum state, at some particular time t = to, can be identified with the state 
that defines the effective potential through the static variational principle 
(Eqs. (Hl-OBD). 



4 An application. The harmonic oscillator 

We now discuss a simple application of the results presented in the previous 
section. As in [15], we consider a trivial field theory (in one time and zero 
space dimension), namely the case of a one-dimensional harmonic oscillator 
whose Hamiltonian is {% = m = 1): 

A= p 

2 2 2 K J 

The term — ^ has been added in order to have a vanishing vacuum energy. 

In this case it is known that the effective action T coincides with the 
classical action S. In [15] this result has been recovered by means of the 
variational principle, which we illustrated in Section 2, and the explicit so- 
lutions of Eqs. ((TH) and (fTB*)) for a generic source function J(t) have been 
shown to be: 



(Q\i> + ,t) = (Qty_,t) = (^) V4 exp {-|(Q - qf - iq{Q - q) - *-f } , 

(25) 

with 

q(t) = (VM|QhM) = - / dt' J{t') e— (26) 

Also, the Lagrange multiplier w(t) turns out to be: w(t) = J(t)q(t). 

From Eq. ()25|) we see that \ifi+,t) = \ij)-,t) as was expected because of 
the periodicity of the motion of q for the harmonic oscillator effective action. 
Moreover, we see that the right-hand sides of Eqs. ()25|) and (|2b[) are totally 
determined once the source J is given. 

In addition we show a concrete realization of the procedure illustrated in 
section 3. In fact, we derive from the effective action the equation of motion 
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for q(t) which is equal to the classical one (being in this case S[q] = T[q}) 
and we therefore have the typical solution 



q(t) = A cos(o;t + ur), (27) 

where the constant r is fixed by the boundary conditions. By taking 

nix 

r = - (28) 

with n integer and to = ~ T (where to is defined in Section 3), we get q(t ) = A 
and q(to) = 0. Moreover q(t) covers n periods in the time range [— r, r] and 
so q{r) = A and q{j) = 0. 

Our problem here is to show that the source J(t) that drives q(t) from 
q(— oo) = (in fact as seen in Section 3, the quantity q at t = ±oo must be 
the expectation value of the coordinate for the ground state of the harmonic 
oscillator, which is zero) to the value q{—r) = A. Then, after the interval 
— r < t < t where J = and q(t) has the mentioned periodic behaviour, q 
must evolve from q{r) = A to q(+oo) = 0. 

This source is 

J(t') = A u e~ elt ' 1 (-00 < t' < -r) and (r < t' < +oo) (29) 

J(t') = (-r < t' < r) (30) 

where e is a small parameter, which eventually has to be sent to zero. In 
fact by inserting this particular source into Eq. (j2HJ), q(t), for — r < t < r is 
given by 

q (t) = li m l im / \' T dt' e-^ e^' e e4 ' + T dt' e iuJt e'^' e~ et ' ) 



2 T^+oo 

= A cos(wt + wr), (31) 

where we have introduced the intermediate step of cutting the time integrals 
at a large time T > and then taking the limit T — * +oo, which will be 
helpful in the following calculations. The final result in Eq. (|31|). which is 
the desired solution introduced in Eq. (|2*7j). has been obtained with the help 
of 

exp(— iur) = exp(za;r), (32) 
which holds because of Eq. (J2S|) . 
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We also check the behaviour of q(t) for t — > +00, as obtained from the 
source defined in Eqs. (J29|) and (|3U|) . We therefore compute the integral 
in Eq. ()26|) by introducing as before the cutoff ±T for large (positive and 
negative) values of t' and putting t = T, and then taking the limits T — > +00 
and e — > 0: 



Instead of taking the limit T — > 00 in the last line of Eq. (|33|). we keep T 
fixed to a very large value due to the oscillating terms exp(— itoT) and perform 
the limit e — > 0. Then, by making use of Eq. (|32|) . we find q(+oo) = 0. The 
computation for q(— 00) is totally analogous and gives again g(— 00) = 0. 

Equations (}3*T]) and (J3l!j) provide the desired results for the source J(t) 
defined in Eqs. f!29|) and (j3*Uj) . As noted in Section 3, there is a discontinuity 
in the function J(t) at the two points t = ±r; however, in this particular 
example there is no consequence on the coordinate expectation value (and 
its derivatives), which are continuous at those points. 

According to Section 3 we must also find a new source that determines 
a static coordinate expectation value, namely q(t) = A for — r < t < r. 
For this purpose, we retain the definition in Eq. ()29|) for large positive and 
negative times but, following the argument in Section 3, instead of Eq. (|3T?|). 
we now take J = A to exp(— er) constant in the time interval [— t, t] and 
this new source does not have any discontinuity points. Moreover in the two 
asymptotic regions the integrals for t' are equal to the ones performed above 
and we have only to compute the contributions in the region where the former 
source vanishes according to Eq. (j3TI|) . Therefore, for — r < t < r, q(t) is 
given by Eq. (J31j) plus an additional term (note that, when performing an 
integral in this limited range [— r, r] the term exp(— er) in the new source will 
always tend to 1 for e — > and it can therefore, for simplicity, be neglected 
from the beginning): 





(34) 
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where again Eq. has been used. This is exactly the required solution 
for q. 

Finally, in order to check that q(t) vanishes asymptotically even with this 
source, we must check that the vanishing result in Eq. ()33j) is not modified. 
This time we have the additional contribution 

q{+oo) = J T dt > e -MT~t') = _ e -i^ = Q; (35) 

which is vanishing because of Eq. (|32j) . The same is valid for q{— oo). 

5 The double-well potential 

As a further test of the arguments presented in Section 3, we now consider 
the non-trivial case of the motion of a wave packet in a double-well potential. 
To follow the dynamical evolution of the position operator expectation value, 
we make use of Eq. (Q) within the framework of the LPA (see Eq.fjHJ)) for the 
effective action. The effective potential V e ff(q) is obtained with the help of 
the wilsonian RG equation (see below), which, as we mentioned before, gives 
a non-perturbative, convex, effective potential (for comparison also the one- 
loop effective potential will be considered at the end of this section). Again 
we fix for q and q the two initial conditions: 

q(to) = qo and q(t ) = 0. (36) 

Then, being again in the presence of a periodic motion, we are allowed to 
use Eq. (JTJ) as an equation of motion for q(t) = (Q) t = (ip,t\Q\ip,t). 

According to our arguments, the function q(t) derived by means of the 
above procedure should be equal to the one obtained by solving the time- 
dependent Schrodinger equation with the form of the initial wave packet fixed 
by Eq. (j22j) . In fact we have numerically solved the latter problem with the 
initial conditions given by Eq. (}2*2*j) : for comparison, we also considered 
the evolution of an initial gaussian wave packet that is often used in the 
applications 4 . 

When we compare these results with the one obtained within the frame- 
work of the effective action formalism, we see that, given the initial conditions 

4 The RG flow equation for Vrg and time-dependent Schrodinger equation were solved 
with the help of the NAG routines [27]. These routines were also used to compute the 
mean values of Q and P. 
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Figure 1: V dw (q) and V RG (q) for A = 0.15. 



(jHHJ), Eq. ((H) (together with Eq. (JHJ)) provides a very accurate description of 
the dynamical evolution of q(t), if the initial state is chosen according to our 
criterion, i.e. with the help of Eq. (J22|)- We also see that the (non-convex) 
one-loop effective potential, that is the approximation to V e ff considered in 
the applications, is inadequate to describe the dynamics of q(t). 

In Fig. [U we have plotted the double-well potential Vd w (q), which we have 
written as: Vd w (q) = — \q 2 + Xq A + together with the effective potential 
V e ff(q), obtained by solving the RG equation [28-31]: 



(here U k (q) means g p q ' ) . Note that in this framework the classical double- 
well potential Vd w (q) is nothing but the bare potential, i.e. the UV boundary 
condition for Uk(q); also note that the effective potential V e ff(q) is approx- 
imated by the solution of Eq. (|3*7|). once it has been integrated down to 
k = 0, i.e. once the quantum fluctuations have been taken into account: 
Veff{o) — ^rg(Q') = Uk=o(q)- Also, when Eq. Q is inserted in Eq. ((TJ), the 
quantum equation of motion for q(t) takes the form: 




(37) 
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(b) 



4 t 5 



Figure 2: The case A = 1.1. (a). The initial wave packets ip ,{q) and ip (q)- 
(b). The time evolution of q(t) as obtained from the Schrodinger equation 
with initial wave packets ipj(q) and ip g {q), and from Eq. ()38|) with the initial 
conditions given in the text. 



Q = -fijV 9ff (q) , (38) 

i.e. the dynamical evolution of q(t) is given by the classical equation of 
motion, being the classical potential replaced by V e ff, which, as discussed 
above, will be approximated by Vrq m our numerical analysis. 

We now illustrate the above points by computing q(t) and p(t) (= q(t)) 
for three different values of A : 1.1 , 0.15 and 0.07. We also choose the initial 
time as to = 0, and the initial values of q and q as follows: g(0) = 0.2 for 
A = 1.1, g(0) = 0.5 for A = 0.15, g(0) = 1.0 for A = 0.07, and g(0) = for 
each value of A. 

According to our arguments, for such a choice of the initial values, Eq. (JTJ) 
describes the time evolution of the expectation value q(t) related to the initial 
wave packet, which is obtained by solving the time- independent Schrodinger 
equation (see Eq. (J'22|)): 
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Figure 3: The phase- space q — p diagram for the three cases considered in 
Fig. 2b (A = 1.1). 

-—A + V dw {q) -Jqj iP(q) = E J ^{q) , (39) 

where, for each of the chosen initial values of q(t), J = J(q ) is computed 
with the help of the equation: (il> j\Q\ijj j) = q (for notational convenience 
we have replaced g(0) by q ). 

Once the initial wave packet ip(q, 0) is obtained with the help of Eq. 
dHSJ), ip(x,t) is obtained by solving the time-dependent Schrodinger equation 
numerically, and the exact values of q(t) and p(t) are then computed. These 
results are compared with those obtained by solving Eq. ((TJ). 

As we already discussed in the Introduction, the boundary conditions 
that define the effective action are not appropriately taken into account in 
the applications of the effective action formalism to dynamical problems . 
Typically one associates a gaussian initial wave packet with the initial con- 
ditions (JHEI)- For this reason we have also considered the exact evolution of 
a gaussian wave packet, ip g (q), chosen in such a way that its width is equal 
to that of the ground state of the harmonic approximation at the bottom of 
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Figure 4: The case A = 0.15. (a). The initial wave packets ipj(q) and ip g (q). 
(b). The time evolution of q(t) as obtained from the Schrodinger equation 
with initial wave packets ipj(q) and ifi g (q), and from Eq. (JHHj) . with the initial 
conditions given in the text. 



one of the two wells. 

For each of the three values of A considered, we present in the following 
figures the shape of the initial wave packets, ipj{q) and ip g {q), the expectation 
value of the position operator (Q) t = q(t) as a function of time, and the 
phase-space, q — p, diagram. 

In Figs. El and El we present the results for the A = 1.1 case. From Figs. 
Eb and El we see that the results obtained with tpj(q) and 4> (q) are almost 
the same (and that they are both well approximated by the results of Eq. 
(|38|) V As is clear from Fig. |2t, this is because, for this particular value of A, 
il) g {q) almost coincides with il>j(q)- 

However, we shall now see that, when we consider other values of A, 
corresponding to different heights of the potential barrier, the choice of the 
correct initial wave packet becomes crucial. Taking for instance A = 0.15, we 
obtain two different shapes for if)j(q) and ifi g (q) (Fig. 0^,). It is not difficult 
to imagine that they have a different dynamical evolution, as can actually 
be seen from Figs. EJd and 03 The crucial result for our analysis is that Eq. 
()38|). i.e. the LPA of Eq. gives a very good approximation to the time 
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Figure 5: The phase-space diagram for the three cases considered in Fig. 4b 
(A = 0.15). 

evolution of (Q) when the initial state is tpj(q) (see Fig. EJ). This gives a 
very robust support to the arguments we have developed in Section 3. 

As an additional example, we have considered the case A = 0.07. Once 
more, Figs. El and [7| confirm that the effective action formalism, in particular 
the LPA we have considered in the present paper, actually describes the 
dynamical evolution of the expectation value of the position operator once 
the initial state is selected according to Eq. 

It is worth noting that, because of the initial conditions chosen for q, 
the cases considered above describe tunnelling processes. Needless to say, 
a semiclassical expansion of T at the lowest order, i.e. the approximation 
T — S, would have been unable to describe these tunnelling processes. In 
fact, in this approximation, V e ff — Vd w . 
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Figure 6: The case A = 0.07. (a). The initial wave packets ipj(q) and ip g (q). 
(b). The time evolution of q(t) as obtained from the Schrodinger equation 
with initial wave packets ipj(q) and ip g {q), and from Eq. ()38|) with the initial 
conditions given in the text. 



From the above results, it is clear that the trajectory of the gaussian wave 
deviates more and more from the correct result for decreasing values of A, 
i.e. for larger quantum effects due to the tunnelling. Therefore it can be 
employed as an ansatz of the correct wave packet only for a specific range 
of the coupling A and this is exactly what has been done in [5], where the 
results obtained with a gaussian wave packet are still accurate (as is the case 
for A = 1.1 as considered above). 

This is a clear indication that a semiclassical approximation to Eq. 
which, for perturbatively small values of A, would correspond to the periodic 
motion of an almost gaussian wave packet around one of the classical minima, 
is misleading. The trajectory obtained from Eq. ([TJ includes quantum effects 
that are neglected in a semiclassical approach. 

Figure [7| needs one more comment. In fact for this smaller value of A the 
phase-space q — p diagram obtained with the "correct" initial wave packet, 
even though it is still qualitatively well described by the phase-space diagram 
obtained with Eq. (J38|) . shows a certain deviation from the latter. The reason 
for such a deviation has already been explained in [5] and has to be traced 
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Figure 7: The phase-space diagram for the three cases considered in Fig. 6.b 
(A = 0.07). 

back to the approximation that we have used for T, the LPA. Actually, it 
was shown in [5] that the successive step in the derivative expansion of the 
effective action, i.e. the inclusion of the wave-function renormalization in 
Eq.(jnj), provides an improvement of the LPA results. Here we are concerned 
with different problems, and so no longer pursue this issue. 

Before ending this section we would like to add a comment concerning 
the approximation of V e ff, namely the one- loop approximation: V e ff ~ Vn. 
In Fig. |HK we show the one-loop effective potential versus the RG effective 
potential. As is well known, the one-loop potential does not enjoy the con- 
vexity property of the exact effective potential K//- Actually this is true at 
any finite order of the loop expansion for the effective potential. Moreover, 
in the region between the two inflection points of the classical potential, the 
one-loop potential develops an imaginary part; according to the usual inter- 
pretation [21], this signals the instability of the configurations in this region. 
In Fig|H^L the real part of the potential is shown, which is the only part that 
is taken into account in dynamical problems. 
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Figure 8: (a). The one-loop (non convex) effective potential compared with 
the RG (convex) one. (b). The phase-space diagrams obtained from Eq. 
(J58J) with V e ff approximated by Vu and Vrg, for a motion with q(Q) = 0.5 
corresponding to an energy below the potential barrier. Differently from Vrg, 
Vu completely fails in describing the tunnelling. In (a) and (b) it is A = 0.15. 



As is well known, Vu can be obtained analytically by replacing U k by U K 
in Eq. ()37|) (remember that U\ = Vd w ), and then performing the elementary 
integration from k = to k = A, i.e. by solving the RG flow equation in 
the independent-mode approximation. We thus obtain: Vu = Vd w + SVu. 
Expanding in 1 and neglecting ^-independent terms, we can write the one- 
loop correction 5Vu as : 



SVu = yvZ arctan^= + 0(1) for|x| > 

(40) 



'12A 



From the above equations we see that, in the region \x\ < 77^, T^(SVu) 
vanishes with A — > 00. A comparison of Fig. |H^, with Fig. ^also shows that, 
within this inner region, Vu is practically equal to Vd w . This is an example 
of the well-known fact that the loop expansion cannot change the concavity 
of Vdw 
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Figure 9: The phase-space diagrams obtained from Eq. (|H8j) . with V e ff 
approximated respectively by Vu and Vrg, for a motion with q(0) = 0.7 
corresponding to an energy above the potential barrier. Even if now also Vu 
allows the whole g-region to be explored, again we see that the motion is 
quite well described by Vrg an d not by Vu. Here we have taken A = 0.25. 



It is clear that for the kind of problems we have considered in this section, 
namely for tunnelling processes, the one-loop effective potential Vu is as 
inadequate as the classical potential Vd w to approximate V e ff in Eq. (jSHJ). 
This is illustrated in Fig. \£jp, where it is shown that for an initial energy 
below the energy barrier the motion is obviously confined within one well. 
According to the results of the semiclassical expansion, we can actually state 
that at no finite order of this expansion can these phenomena be described 
by replacing the semiclassical effective potential in Eq. (|38|) . Instead they 
are perfectly well described by the solution of Eq. (JH7j) . Vrg, which is only 
the lowest-order approximation to V e ff of a different expansion, namely the 
gradient expansion. 
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In order to push this comparison a step further, we have also considered 
a motion with energy above the potential barrier. Naturally, in this case, by 
replacing Vu in Eq. (JHHj) we obtain a solution that, as is the case for Vrg, 
explores the whole allowed region in position space. However, as is shown in 
Fig. El the correct description of the motion is again given by Vrg and not 
by Vu- 

6 Summary and outlook 

In the present work we have shown that in order to make use of the equation 
3^7]^ = (Eq. (JTJ) in the text) as a "quantum equation of motion", attention 
has to be paid to a certain number of important issues. 

First of all we need to make sure that we are under conditions such that 
the argument of T, the classical field (fi(x, t), is actually the expectation value 
of the quantum operator <3>(x) in a given state t), which is not always the 
case. 

Once we are under these conditions, we have shown that the asymptotic 
|in) and |out) states that enter the definition of the effective action have to 
be taken appropriately into account in the determination of the initial wave 
functional ^[(J)(x, to)} associated with the motion of (p(x, t). More precisely, by 
considering a derivative expansion for T, we have argued that it is possible to 
determine uniquely the initial state ^>[(f)(x, to)} from the boundary conditions 
for the motion of (p(x,t), assigned for instance on a given manifold t = to 
(initial time), together with the asymptotic conditions encoded in the |in) 
and | out) states. 

For initial conditions that are particularly important in the applications, 
namely the case of constant initial <fi and vanishing derivatives of <fi at t = t Q , 
we have been able to show that the initial (functional) wave packet obeys a 
modified time-independent Schrddinger equation, actually the equation that 
allows the definition of the effective potential. 

Finally we have shown that, in the framework of the LPA for T, where 
the quantum equation of motion is obtained by the classical equation of 
motion once the classical potential is replaced by the effective potential, 
a reliable approximation to the effective potential comes from the solution 
of the corresponding wilsonian RG equation; instead, the one-loop effective 
potential, and more generally any finite-order approximation to the loop 
expansion for K//, fails to describe the dynamics of (f>(x,t). 
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All the above points have been illustrated by considering quantum me- 
chanical examples. After a brief application of these results to the simple 
case of the harmonic oscillator, which confirmed our arguments, we have 
considered the case of a double-well potential. The nice feature of consider- 
ing a quantum-mechanical example is that we have the possibility to compute 
also the dynamical evolution of the wave packets with the time-dependent 
Schrodinger equation, i.e. we can compare our results with exact results. 
The cases we have investigated all confirm the correctness of our conjectures. 

One point that was not considered in the present work concerns the pos- 
sible extension of these results to the case of higher orders in the derivative 
expansion of T. It would also be interesting to investigate the impact of the 
presence of non-local terms in T. We plan to investigate this point in future 
work. 

Naturally our ultimate goal is to apply our results in a QFT context 
where, as we cannot solve the quantum mechanical (Heisenberg or Schrodinger) 
equations exactly, we really need to have sensible approximation schemes. We 
hope that our results can help in this direction. 



26 



References 

[1] G. Jona-Lasinio, Nuovo Cimento 34 (1964) 1790. 

[2] R.J. Rivers, Path Integral Methods in Quantum Field Theory, Cam- 
bridge University Press (1987). 

[3] Lowell S. Brown, Quantum Field Theory, Cambridge University Press 
(1995). 

[4] See for instance: J. Zinn- Justin, Quantum Field Theory of Critical Phe- 
nomena, Oxford Science Publications (1993). 

[5] G. Andronico, V. Branchina and D. Zappala, Phys. Rev. Lett. 88 (2002) 
178902. 

[6] R. Ramirez and T. Lopez-Ciudad and J.C. Noya, Phys. Rev. Lett. 81 
(1998) 3303. 

[7] R.H. Brandenberger, Rev. Mod. Phys. 57 (1985) 1. 

[8] D. Boyanovsky, H.J. De Vega, R. Holman Phys. Rev. D49 (1994) 2769. 

[9] F. Cooper, S.-Y. Pi and P. Standoff, Phys. Rev. D34 (1986) 3831. 

[10] S.-Y. Pi and M. Samiullah, Phys. Rev. D36 (1987) 3128. 

[11] B. De Witt, The Global Approach to Quantum Field Theory, Oxford 
University Press (2003). 

[12] J. Schwinger, J. Mat. Phys. 2 (1961) 407. 

[13] P. M. Bakshi and K. T. Mahantappa, J. Mat. Phys. 4 (1963) 1. 

[14] E.Calzetta and B.L.Hu, Phys. Rev. D15 (1987) 495. 

[15] R. Jackiw and A. Kerman, Phys. Lett. A71 (1979) 158. 

[16] S. Coleman, Aspects of Symmetry, Cambridge University Press (1985), 
pp. 132-142. 

[17] K. Symanzik, Commun. Math. Phys. 16 (1970) 48. 

[18] T. Curtright and C. Thorn, J. Math. Phys. 25 (1984) 541. 

27 



[19] D.J. Callaway and D.J. Maloof, Phys. Rev. D27 (1983) 406; D.J. Call- 
away, Phys. Rev. D27 (1983) 2974. 

[20] V. Branchina, P. Castorina and D. Zappala, Phys. Rev. D41 (1990) 
1948. 

[21] E. J. Weinberg and A. Wu, Phys. Rev. D36 (1987) 2474. 

[22] R. Fukuda, Prog. Theor. Phys. 56 (1976) 258. 

[23] A. Ringwald and C. Wetterich, Nucl. Phys. B334 (1990) 506. 

[24] N. Tetradis and C. Wetterich, Nucl. Phys. B383 (1992) 197. 

[25] J. Alexandre, V. Branchina and J. Polonyi, Phys. Lett. B445 (1999) 
351. 

[26] A. Horikoshi, K.-I. Aoki, M. Taniguchi and H. Terao, |hep-th/9812"050l 
in Proc. Workshop on the Exact Renormalization Group, Faro, Portu- 
gal, 1998, World Scientific, Singapore (1999); A.S. Kapoyannis and N. 
Tetradis, Phys. Lett. A276 (2000) 225; D. Zappala, Phys. Lett. A290 
(2001) 35. 

[27] The Numerical Algorithms Group Limited (1997), NAG Fortran Library 
Manual, Mark 18. 

[28] F. Wegner and A. Houghton, Phys. Rev. A8 (1973) 401. 

[29] A. Hasenfratz and P. Hasenfratz, Nucl. Phys. B270 (1986) 687. 

[30] A. Bonanno and D. Zappala, Phys. Rev. D57 (1998) 7383. 

[31] A. Bonanno, V. Branchina, H. Mohrbach and D. Zappala, Phys. Rev. 
D60 (1999) 065009; V. Branchina, Phys. Rev. D62 (2000) 065010. 



28 



